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Abstract. Small-angle correlations of pairs of nonidentical light charged particles produced in central 
■ collisions of heavy ions in the A = 100 mass region at a beam energy of 400 yl-MeV are investigated with 

$H ' the FOPI detector system at GSI Darmstadt. The difference of longitudinal correlation functions with the 

relative velocity parallel and anti-parallel to the center-of-mass velocity of the pair in the central source 
frame is studied. This method allows extracting the apparent space-time differences of the emission of 
the charged particles. Comparing the correlations with results of a final-state interaction model delivers 
quantitative estimates of these asymmetries. Time delays as short as 1 fm/c or - alternatively - source 
radius differences of a few tenth fm are resolved. The strong collective expansion of the participant zone 
introduces not only an apparent reduction of the source radius but also a modification of the emission 
times. After correcting for both effects a complete sequence of the space-time emission of p, d, t, ^He, a 
particles is presented for the first time. 
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1 Introduction 

Correlations of two nonidentical particles at small relative 
velocities are - due to final state interactions - sensitive 
to the space-time structure at freeze-out. Usually, the cor- 
relation is used to extract the size of the source and the 
time duration of the emission However, it 

contains also information on the emission time differences 
of the two particles. Gelderloos and Alexander proposed 
to construct velocity difference spectra at small relative 
angles ||]. Comparing these spectra with results of trajec- 



tory model calculations they were able to infer the emis- 
sion order and the time intervals between the emission of 
the two particles If the faster particle of a pair ap- 
proaches and passes the other one from behind the pair 
experiences a stronger final-state interaction than in the 
case that the faster particle has started in front of the 
slower one. Thus, the ratio of the correlation functions 
with the relative velocity parallel and anti-parallel to the 
pair velocity is related to the formation sequence. 
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For heavy-ion collisions of a few hundred MeV per nu- 
cleon mainly nucleons and light composite particles with 
charge number Z < 2 are emitted, and one expects rel- 
ative short time differences. For these light ejectiles clas- 
sical trajectory calculations are not accurate as quantum 
effects dominate at momenta larger than Ti/aQ where ao is 
the Bohr radius. Very recently, Lednicky et al. have 
proposed to use the above sensitivity of the correlation 
function to directly measure the space-time differences in 
the emission of particles of different types. This is possi- 
ble since the wave functions for nonidentical particles are 
asymmetric with respect to the forward and backward di- 
rection of the relative momentum. First theoretical |0 
and experimental work |l2) concentrated on the determi- 
nation of the asymmetry of particle production of proton- 
pion pairs. In the present work the proposed method is 
applied to measure the pair-wise space-time differences of 
different light charged particles. Applying directional cuts 
on all relative- velocity correlation functions of nonidenti- 
cal particles allows one to determine - even with a certain 
redundancy - the whole sequence of space-time emission 
points of p, d, t, '^He, and a particles. 

The two-particle correlations are sensitive to the spa- 
tial separation of the two particles at the time when the 
second of them freezes out. This separation is the sum of 
the spatial separation between the two freeze-out points, 
and the time separation, multiplied by the velocity of the 
first particle. A correlation measurement does not allow to 
disentangle these two components. Therefore, in the fol- 
lowing investigation we will discuss two cases, first that 
all particles are emitted from the same region and second 
that they emerge at the same time instant. 



2 The experiment 
2.1 Detector setup 

The experiment has been performed at the heavy-ion syn- 
chrotron SIS at GSI Darmstadt. Targets of 1 % interaction 
thickness of ^^Ru and ^^Zr have been irradiated by ^^Ru 
and ^^Zr ions of 400^-MeV beam energy. In order to get 
sufficient statistics, the data of all target-projectile combi- 
nations have been used. The original aim of the utilization 
of target and projectile nuclei with equal mass but differ- 
ent isospin content was to answer the question whether 
the colliding nuclear system attains a full thermo-chemical 
equilibrium during the collision process. The first exper- 
imental results reveal substantial transparency effects in 
phase space regions already slightly apart of midrapidity 

The present analysis uses a subsamplc of the data, 
taken with the outer Plastic Wall/Helitron combination 
of the FOPI detector system |l4|. The Plastic Wall de- 
livers - via energy loss vs. time-of-flight (TOF) mea- 
surement - the nuclear charge Z and the velocity (3 of 
the particles. The Helitron gives the curvature (which 
is a measure of the momentum over charge {p/Z)) of 
the particle track in the field of a large superconducting 
solenoid. Since the momentum resolution of the Helitron is 
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Fig. 1. Two-dimensional distribution of yields (fa/dptdy of 
p, d, t, a particles in the pi — if plane for central reactions 
selected by a 8 % cut on large charged particle multiplicities. 
Target and projectile rapidities are given by = -1 and -1-1, 
respectively. The full lines are levels of constant yield of 20, 40, 
60, and 80 % of the maximum value. Dashed lines represent 
the polar angle limits at 8.5 and 26.5 degrees. 

rather moderate, this detector component serves for par- 
ticle identification only. The mass m is determined via 
mc = {p/Z)Heil{Pl/Z)piaWa. where 7 = (1 " Z?')"'/". 
The Plastic Wall and the Helitron have full overlap only 
for polar angles between 8.5 degrees and 26.5 degrees. 
The corresponding flight paths amount to 450 cm and 
380 cm, respectively. Monte-Carlo simulations have been 
performed in order to study the influence of the finite de- 
tector granularity and of the TOF and position resolutions 
on the velocity and finally on the proton momentum. The 
resolution of both quantities is governed by the TOF reso- 
lution, which is axoF = 80 (120) ps for short (long) scin- 
tillator strips located at small (large) polar angles [^ . 
The detector granularity delivers a negligible contribution 
to the velocity resolution jl5|. Thus, between midrapid- 
ity (j/cm = 0.447, Pern = 0.419) and projectile rapidity 
{Vproj — 0.894, (3proj — 0.713) the velocity can be deter- 
mined with a precision of (Tp/ (3 ~ 0.4% — 0.8%. Finally, 
from the velocity other kinetic quantities like the velocity 
^cm ^j^g particle momentum p"^™ after transformation 
into the cm. system of the colliding nuclei are deduced. 



2.2 Event classification 

About 4-10^ central collisions are selected by demand- 
ing large charged-particle multiplicities to be measured 
in the outer Plastic Wall. The corresponding integrated 
cross section comprises about 8 % of the total cross 
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section. For this centrality class one would expect - 
within a geometrical picture - an average impact pa- 
rameter of about 2 fm. Simulations which we have per- 
formed with the IQMD model jl^ predict an average 
impact parameter of about 2.5 fm. For different par- 
ticles with mass number Adus Fig.^ shows the phase 
space coverage of the detector components, outer Plas- 
tic Wall/Helitron, in the transverse momentum vs. rapid- 
ity plane for events selected by this centrality condition. 
Here, p° = {pti A.i^^) I [Ppro, / Aproj)cm = {Ptl)l{l3lYj:::' 
and ?/° = {y/ypro])cm = {y/Vcm - 1) are the normalized 
transverse momentum and rapidity, respectively. Both 
observables are related to the corresponding projectile 
quantities in the cm. frame of the colliding nuclei (with 
{PlYcm' = 0-462 for 400 A-MeV beam energy). It is ob- 
vious that for central collisions the Plastic Wall preferen- 
tially measures midrapidity particles with small velocities 
in the cm. system {{v""^) ~ 0.25 c - 0.30 c). 

In previous investigations of central Au-|-Au collisions 
between 100 and 400 A-MeV beam energy it was found 
that the correlation function of pairs of intermediate mass 
fragments (IMF) is strongly affected by the collective di- 
rected sideward flow of nuclear matter ||r^,|I^. This di- 
rected sideflow causes an enhancement of correlations at 
small relative momenta. The enhancement results from 
mixing of differently azimuthally oriented events; it van- 
ishes if the events are rotated into a unique reaction plane, 
which is determined by the standard transverse momen- 
tum analysis Q. This procedure is allowed because the 
geometrical acceptance and the detector efhciency are az- 
imuthally symmetric. Thus, the technique of event rota- 
tion is applied also to the present data in order to prevent 
that such artificial correlations are introduced into the ref- 
erence momentum distribution of the correlation function 



(cf. Sect. 2.3) 



2.3 Correlation function 

Let Yi2(pi,p2) be the coincidence yield of pairs of par- 
ticles having momenta pi and p2. Then the two-particle 
correlation function is defined as 

1 , T-i / \ kt ^-^events, pairs ^2 (pi , P2) . 

l-t-R(pi,P2) =AA= 7——Y (1) 

2^ events, pairs 12,mixWl, V2 ) 

The sum runs over all events fulfilling the above mentioned 
global selection criterion and over all pairs satisfying cer- 
tain conditions given below. Event mixing, denoted by the 
subscript " mix" , means to take particle 1 and particle 2 
from different events. A/" is a normalization factor fixed 
by the requirement to have the same number of true and 
mixed pairs. The statistical errors of all the correlation 
functions presented below are governed by those of the 
coincidence yield, since the mixed yield is generated with 
two orders of magnitude higher statistics. The correlation 
function (1) is then projected onto the relative momentum 

q, 

q = MVi2-M(vr-vr)- (2) 



Here, v^™ are the particle velocities calculated in the cm. 
system of the colliding nuclei and — [17111712) / {mi + m2) 
is the reduced mass of the pair. Besides the above de- 
scribed global event characteristics we use gate conditions 
on the angle C between q and the cm. sum momentum of 
the particle pair Pf™ = pf™-|-p2™ and on the pair velocity 
V = |Pf™|/(mi-f TO2). In order to exploit the full available 
statistics, two complement types of longitudinal correla- 
tion functions are generated. The forward and backward 
correlation functions are defined by cuts on the angle Q, 
cosC > and cosC < 0, respectively. This choice selects 
pairs with the longitudinal velocity component vl (pro- 
jection onto the pair velocity) of particle 1 being greater 
or smaller than the corresponding value of particle 2: 

R+(g) = l + R(g,cosC>0) = l + R((z,VL,i > VLa) (3) 

R- {q) = l + R(9, cos C < 0) = 1 + R(9, VL,i < VL,2) (4) 

From the velocity resolution as estimated in Sect.^ the 
corresponding q resolution is deduced. It is expected to 
amount to 6q = /i(5wi2, where 6vi2 = ■\/2(((5w'^™)^) = 
(0.006 ± 0.002) c 

3 Analysis 

3.1 Correlations from final-state interaction 

The correlation function of tw o par ticles 1 and 2 which 
move with a pair velocity V is [|3|,|20| 

l + i?(pi,P2) = J dtidt2dridr2Pi{'V,ri,ti)p2{V,r2,t2) x 

|>f'q(ri-r2-V(ti-i2))P, (5) 

where the density /9i,2(r, i) describes the probability to 
find particle 1 or 2, respectively, at time t and spatial co- 
ordinate r from which they are emitted (freeze-out config- 
uration) . The wave function !Z'q describes the relative mo- 
tion of the two particles with momenta pi and P2. Since 
the interaction with the source is neglected this wave func- 
tion depends only on the relative coordinate r = ri — r2 
and on the relative momentum q defined in (^. Here, 
we assume that the source functions pi.2 do not depend 
on the velocity and have a Gaussian shape in space and 
time characterized by radius parameters Ri,2 and emis- 
sion times ri.2 while the duration of the emission is given 
by To: 



Pi,2(r,i) 



iTT^l 2^0 



exp[ 



2tI 



(6) 



The simple form of Eq. (g) allows us to integrate over the 
center-of-mass coordinate (in\T\ +m2V2)/{mi -I-7712) and 
the time variables leading to 



i?(pi,P2) = J dV^(r,V)|tf'q(r) 



(7) 
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with 



5(r,V) = 



(47r)-V2i?2£, 



cxp 



1 



{YAt + r 
' (r-Vf) 



(8) 



where we have introduced the time difference Ati2 and an 
effective radius D by 

At,2 = n-r2, = RIHVto)\ Rl = \{Rl+Rl)- 

(9) 

For a finite emission time difference Ati2 the source 
function 5 depends on the sign of the scalar product V • r. 
This asymmetry with respect to r is transferred to the 
correlation function if the wave function tf'q contains 
terms of an odd power of the scalar product q • r which is 
possible for nonidentical particles. 

The wave function is generated by partial wave 
expansion technique. To incorporate the spin degrees of 
freedom we use either spin-spin coupling (for p-p, d-p, 
t-p, •^He-p, d-d, t-d, ^He-d, and ■^He-t correlations) or 1-s 
coupling (for a-p, a-d, a-t, and a-'^He correlations). Thus, 
the partial waves are classified by angular momentum and 
either total spin s = si -t- S2 or j = / -f S2- There is also 
an option to include partial waves with total angular mo- 
mentum J = / -|- s to describe dominant resonances. 

The corresponding radial Schrodinger equations are 
solved using the Coulomb potential and Woods-Saxon po- 
tentials 



V{r) 



1 -I- exp[-(r - Rws)/an,s] ' 



(10) 



For each partial wave, the parameters Rws^ o,ws, V^js are 
chosen such that the phase shifts are reproduced. As al- 
ready mentioned in ref. ||^ it is important to find poten- 
tials which generate the correct dependence dS/dq since 
the derivative of the phase shift is the relevant quantity 
in determining whether correlations are suppressed or en- 
hanced. In most of the cases we make use of the param- 
eters already obtained in refs. j^,|2l|,^. (Note that for q 
values approaching a certain resonance positions qi the 
phase shift has to amount 90 degrees.) For partial waves 
/ > 3 the nuclear part of the potential is unimportant 
and has been neglected. In the case of p-p correlations the 
wave function has been generated using the Reid soft-core 
potential like in the standard Koonin model 

The wave function 'Fq in Eq. (||) is the projection of the 
many body wave function of the two interacting clusters 
on their relative coordinates. The assumption made above 
that this wave function can be replaced by the partial 
waves could be violated for small distances |ri — r2| where 
many body effects are most important. Consequently, the 
potentials which are obtained by fitting the scattering 
phases might not be the optimum choice to produce the 
correlation function. E.g., the d-p correlation functions 
shown in Fig. ^ are too steeply increasing with relative mo- 
mentum q. The agreement to the measurement is largely 
improved using the potential depth of 6 MeV instead of 



the value of -13 MeV given in ref. § for the s = 3/2, / = 1 
partial wave. However, it is important to note that by 
changing this value the ratio R^/R^ remains essentially 
unchanged. Therefore, we have decided to use potentials 
which are compatible with scattering phases avoiding the 
ambiguity which could arise by searching for new param- 
eters in the very large parameter space. 



3.2 Effect of radial flow 

From Eqs. (||) and (^) we observe that differences of the 
source extensions do not enter in the correlation function. 
This fact is a consequence of the isotropic distribution of 
the emission points which is independent of the particle 
momentum. However, in central nucleus-nucleus collisions 
a considerable amount of the bombarding energy is con- 
verted into fiow energy. This collective expansion of nu- 
clear matter sets in after the compression phase and can be 
observed as a decrease of the slope of the kinetic energy 
spectra. It causes a strong correlation between particle 
momenta and emission points. Therefore, we investigate 
in this section the influence of the flow on both the source 
radius and the emission time. 

For this purpose we introduce into the emission func- 
tion (^) a radial fiow velocity. There is good experience 
with the "nuclear Hubble scenario" which suggests a lin- 
ear velocity profile. Such an assumption has furthermore 
the advantage that the function S{r, V) can be calculated 
analytically when using Gaussian density profiles. Thus, 
we write for particles with different mass numbers Ai^2 



Pi,2(v, r, i) = iVexp[ 



2^1,2 



(t - ri,2f 



^i,2Wo(v - rjvf- 
2T 



(11) 



where TV is a normalization factor which ensures that the 
density is normalized to unity. In ( pi] ) the velocities of the 
particles are thermally distributed, characterized by the 
temperature T, around a radial flow velocity defined by 
the scaling parameter rj. 

The "nuclear Hubble constant" ?/ is connected with the 
radial fiow energy per nucleon 



efio 



moT] 



moT] 



3Rq. 



(12) 



The brackets imply averaging over the Gaussian density 
distribution, and the quantities toq and Rq represent the 
nucleon rest mass and the mean radius, respectively. 

Repeating the integration of Eq. (|^) with the new func- 
tions pi^2 one arrives at a similar relation for the source 
function S as Eq. (||), however the parameters i?o, D and 
Ati2 of Eq. (^) have to be replaced by quantities indicated 
by an asterisk: 



(13) 
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and 



with 



and 



D* = 



_flow 



'1,2 



Tl 2 



1,2 



Al,. 
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A flow _ 

^^12 ^ 



1 + ef{Ai 



1 + £^2^2 



(14) 
(15) 

(16) 
(17) 



Here, we use the abbreviations fi.2 — ^1.2/^Oj a-nd e = 
i flow / itherm is the ratio of the radial flow energy efiow and 
the energy of the random thermal motion ttherm = f T'- 

From Eq. (|l^) one finds that the apparent source ra- 
dius decreases monotonously with increasing energy ratio 
and particle mass. These observations do well compare 
with recent results of the investigation of the sensitivity 
of the proton-proton correlation to collective expansion in 
central Ni-|-Ni collisions at 1.93 A-GeV beam energy p^ . 

The essential result is that even if the two particles are 
produced at the same source region, Ri = R2, and at the 
same time, ti — T2, the flow still causes an effective time 



difference AtT. 



flow 
^12 



if the two particle masses are 



different. This can be understood as follows: For a given 
pair velocity V the heavier partner 1 with the smaller ther- 
mal velocity moves approximately with the flow velocity. 
This means that at freeze-out the heavier particle is lo- 
cated in the region around ri = V/77. Since the position 
of the lighter partner with the larger thermal velocity is 
Gaussian distributed around the center of the source, its 
mean location r2 is behind its heavier partner. Therefore, 
to have a strong final state interaction the lighter needs 
to start earlier. If they start however at the same time it 
seems as if the heavier particle would have started earlier, 
At^2 < in agreement with Eqs. ( |T^ ) and (p^. 

Furthermore, temporal and spatial differences in the 
source distribution add up to the effective emission times 
r* in Eq. (16) quite independently. Therefore, measuring 



the quantity Ati2 alone will not allow us to disentangle 
the two different contributions. In the following analyses 
we will therefore differ two cases, namely first that all 
particles are emitted from the same source size {Ri = i?o) 
and second that they are emitted at the same time (r^ = 0) 
instant. 

In our calculations which take the collective expan- 
sion into account we use an energy ratio of e = 1 and a 
temperature of T = 37 MeV. These values are compatible 
with values obtained in systematic flow studies p3[ of cen- 
tral Au+Au collisions in case of 400 ^-MeV beam energy. 
It should be mentioned that the flow correction zir/j"™ 
changes only slowly with the energy ratio e because the 
temperature T in Eq. (|l^) is correlated with the flow en- 
ergy Cfiow due to the requirement of energy conservation. 
Thus, even if for the present Ru-|-Ru system the ratio e 
would be smaller than for the Au-|-Au system, one would 
find that the correction changes only marginally. E.g. a 
drastic reduction of e by a factor of two leads to typical 
changes of Arf^"^ by about 20% only. 



4 Results and comparison with model 
4.1 The source extension 

Fig.|2| (left panel) shows the angle-integrated two- proton 
correlation function. As in previous fragment-fragment 
correlation analyses [0,|l^ an enhanced coincidence yield 
at very small relative angles is observed, which is due to 
double counting caused mainly by scattering in the scin- 
tillator strips. This disturbing yield is reduced drastically 
by the requirement to match the particle hits on the Plas- 
tic Wall with the corresponding tracks in the forward 
drift chamber Helitron. However, mismatches of tracks 
and scattering processes especially at the wire planes of 
the chamber can give rise to a small amount of double 
counting, too. The remaining left-over of doubly counted 
scattered particles is eliminated by excluding, around a 
given hit, positions within a rectangular segment of az- 
imuthal and polar angle differences |0i — 02| < 4° and 
16*1 — 6*2 1 < 2°. In order to keep the influence of the exclu- 
sion onto the correlation function as small as possible, the 
same procedure is applied to the uncorrelated background. 
However, GEANT simulations |Q have shown that at 
very small relative velocities vu < 0.03 c a small bias of 
the correlation function cannot be excluded [15|. Though 
ratios of correlation functions as our forward/backward 
relation R^/R^ are expected to be rather robust against 
such biases, the corresponding regions in the correlation 
functions are marked as hatched area and are not taken 
into consideration when comparing the experimental data 
with model predictions. 

The experimental p-p correlation function is compared 
to the model for an emission from a source of zero life- 
time and true Gaussian radius Rq which corresponds to 
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+ exp. data 




<v>Tn = 2.2fm daia 



longitudinal 
transverse ^ 




0.05 0.075 0.1 0.125 
q/(GeV/c) 



5 0.075 c 

q/(GeV/c) 



Fig. 2. Left panel: The experimental correlation function of 
proton pairs from central trigger events (dots). The hatched 
area indicates the unreliable region which may be contami- 
nated by doubly counted scattered particles. The full line rep- 
resents the model prediction for an emission from a Gaussian 
source of zero lifetime and radius Ro = 5.7 fm. The dashed 
and dotted lines are the corresponding results for source radii 
which differ from the optimum one by -0.3 fm and +0.3 fm, 
respectively. 

Right panel: The transverse (dots) and longitudinal (squares) 
two-proton correlation functions compared to the model results 
(dashed and dotted lines) for finite lifetime {V}to = 2.2 fm and 
Gaussian radius _Ro ~ 5.4 fm. 
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an r.m.s. radius of Rr 



= = y/3Ra. The theo- 

retical correlation function is folded with an experimen- 
tal resolution function of Gaussia n sh ape with the dis- 
persion ( 5(g) as estimated in Sect. 2^. The apparent re- 
duction (O) of the source radius due to radial expansion 
effects as described in Sect.|3|^ (see also ref. ||l5|) is esti- 
mated to R*o/Ro = l/Vm = I/V2. The best agreement 
of experimental data and model calculations is found for 
i?o = (5.7 ± 0.3) fm. It is obvious that the appearance of 
the correlation peak at g ~ 20 MeV/c allows the deter- 
mination of the source radius with rather high sensitivity. 
This ^He-resonance is the result of the common action of 
the enhancement due to the attractive nucleonic poten- 
tial and the suppression due to both the mutual Coulomb 
repulsion and the antisymmetrization of the wave func- 
tion. The ambiguity of the space and time extents of the 
proton emitting source can be resolved by constructing 
correlation functions with the relative momentum q be- 
ing either perpendicular (transverse correlation function, 
here defined by | cos^l < 0.5) or parallel (longitudinal cor- 
relation function, here defined by | cos C| > 0.5) to the pair 
velocity V. The result is given in the right panel of Fig.||. 
We find a small suppression of the transverse correlations 
with respect to the longitudinal ones which is consistent 
with model predictions for the emission from a source of 
finite lifetime H,^. The best simultaneous fit around the 
correlation peaks (14 MeV/c < g < 42 MeV/c) of both 
the transverse and longitudinal correlation functions de- 
livers Gaussian parameters of the radius (after flow cor- 
rection) and the emission duration of Rq = 5.4 fm and 
{V)to = 2.2 fm, respectively. Obviously, the lifetime effect 
appears unimportant. 

For the determination of the emission time differences 
in the subsequent section the source radius is kept fixed to 
a value of i?i = i?2 = -Ro = 5.7 fm. For simplification the 
duration of the emission tq is set to zero. A finite value 
leads only to a rescaling of the effective radius in Eq. ( p^ . 

4.2 Emission time differences 

Out of the ten different combinations of pairs of non- 
identical particles which can be constructed from p, d, 
t, ^He, and a particles there exist only two correlation 
functions which do not contain contributions of resonance 
decay products. These resonance- free correlation functions 
of deuteron-proton and ■^He-triton pairs are presented in 
Sect. 4.2.1. All other correlation functions contain reso- 



nance contributions of particle-unbound ground states or 
excited states of heavier clusters decaying into pairs of 
light charged particles. In case of the appearance of strong 
and narrow resonances in the q region of interest, the cor- 
responding theoretical correlation functions are corrected 
for the experimental g resolution similarly as for p-p cor- 
relations. 



4.2.1 Final-state interaction without resonances 

Fig. ^ gives the results of d-p correlations. The upper panel 
shows the forward (full dots) and the backward (open 
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Fig. 3. Upper panel; Forward (full dots) and backward (open 
dots) longitudinal experimental correlation functions of d-p 
pairs. The hatched area indicates the unreliable region which 
may be contaminated by doubly counted scattered particles. 
The full and dashed lines give the corresponding model pre- 
dictions with the time delay of table ^. Lower panel: Ratio 
of forward/backward experimental correlation functions (open 
squares) . The full line represents the ratio of the simulated cor- 
relation functions. The dashed and dotted lines give the model 
predictions for times differences deviating by -0.5 fm/c and 
+0.5 fm/c from the optimum one, respectively. 

dots) correlation functions whereas the lower panel rep- 
resents the ratio of both observables. The resonance-free 
correlation functions show a suppression at low relative 
momenta due to final-state Coulomb and nuclear inter- 
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Fig. 4. The same as Fig.H, but for ^He-t correlations. 
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actions. Obviously, the suppression gets stronger in the 
case when deuterons are faster than protons (full dots in 
Fig.||). This indicates that, on average, deuterons are be- 
ing emitted later than protons. Quantitatively, the com- 
parison with the model predictions yields an optimum 
time delay of zir^^ p = "^d ^ ~ fm/c (full line in the 
lower panel). (The dashed and dotted lines in the lower 
panel of Figs.0 and should give an impression how the 
theoretical ratio R^/R^ alters if the time delay is changed 
by -0.5 fm/c and -t-0.5 fm/c, respectively.) However, if one 
does not take into account the time shift due to the radial 
expansion of the participant zone the emission time dif- 
ference would be determined as an apparent value which 
is much smaller (but still positive). Both the true and ap- 
parent time delays are summarized in table|l|. The given 
errors represent the typical band widths of time delay pa- 
rameters around their optimum values which reproduce 
the experimental correlation fimction ratios R^/R^ rea- 
sonably well. 

Fig.^ shows the •^He-t correlation functions and the 
corresponding forward/backward ratio. Since the masses 
of the two species are practically identical, this correlation 
function is the only one where - in case of equal source radii 
- the flow correction (|^ vanishes and, consequently, the 
true and the apparent time delays are identical. Obviously, 
the '^He particles are emitted slightly later than the tritons 
(about 1-3 fm/c, cf. table |l|). 



4.2.2 Final-state interaction with resonances 

If a resonance contribution is dominating in the two- 
particle yield the ratio R^/R" necessarily is forced to 
unity for relative momenta q approaching the resonance 
value Qi. This is due to the fact that for pure two-body de- 
cays both particles are emitted at the same time (and po- 
sition). Indeed, the experimental data show the expected 
behaviour (see arrow positions in Figs. p|-p^. 

Fig. H shows the forward/backward longitudinal corre- 
lation fimctions of t-p pairs together with the correspond- 



Table 1. The apparent and the true values of the emission 
time difference as determined for all combinations of noniden- 
tical charged-particle correlation functions. The source radii 
are fixed to Ri = Rq. 



light-charged 


apparent 


true 


particle 


time delay 


time delay 


combination 






1 - 2 


(fm/c) 


(fm/c) 


d-p 


1.7± 1 


6.5 ±1 


t-p 


-4.2 ±3 


3±3 


^He - p 


4.5 ± 1 


11.7 ± 1 


Q - p 


-2.6 ±3 


6±3 


t - d 


-1.0 ± 1 


1.4 ± 1 


^He - d 


3.0 ± 1 


5.4 ± 1 


Q - d 


-2.0 ± 1 


1.8 ± 1 


^He - t 


1.7± 1 


1.7 ± 1 


a - t 


-2.9 ± 1 


-1.5 ± 1 


^He - Q 


4.4 ± 1 


3.0 ±1 



R" (VL(t)<VL(p)) 
y^^~\l, + (VL(t)>VL(p)) 
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Fig. 5. Upper panel: Forward (full dots) and backward 
(open dots) longitudinal experimental correlation functions of 
t-p pairs. Positions of relevant resonances are marked as ar- 
rows. The hatched area indicates the unreliable region which 
may be contaminated by doubly counted scattered particles. 
The full and dashed lines give the corresponding model pre- 
dictions with the time delay of table ^. Lower panel: Ratio 
of forward/backward experimental correlation functions (open 
squares). The full line represents the ratio of the simulated 
correlations. 

ing ratio. The correlation function exhibits a broad peak 
which contains the contribution of the 1st excited state of 
^He [E* = 20.21 MeV, J'^ = 0+, T = 0.5 MeV, Fp/F = 1, 
qi = 23.6 MeV/c) as well as the 2nd (E* = 21.01 MeV, 
J"" =0-, F = 0.84 MeV, Fp/F = 0.76, q2 = 41.0 MeV/c) 
and the 3rd one {E* = 21.84 MeV, J'^ = 2", T = 
2.01 MeV, Fp/F ^ 0.63, qs = 53.4 MeV/c). The posi- 
tions of the relevant resonances are marked by arrows. In 
the present case only the lowest three excited states of ^He 
with widths F < 2.1 MeV are indicated. Other states in 
the energy region £;* = 23 - 26 MeV (g = 70 - 90 MeV/c) 
are much broader {F > 5 MeV, Fp/F ~ 0.5). Obviously, 
from the lower part of the figure one would conclude a 
negative time delay At^ . However, the radial flow cor- 

rcction ( |17| ) overcompensates this apparent time difference 
and leads to a positive value (cf. table |l|). Due to the pos- 
sible contribution of higher lying resonances of ^He which 
are not taken into account in the model description, a 
rather large error is appended to the time delay. 

Fig. |6| represents the results for ■^He-p correlations. One 
broad maximum shows up which is related to the particle- 
unbound ground state of '^Li {J^ = 2~, F = 6 MeV, 
go = 66.0 MeV/c). An unquestionably positive time delay 
Z\r3jjg p is deduced from the correlation function ratio 

R^/R~. Here, about 60 % of the true time difference given 
in tablej^ arise from the radial flow correction (p^). 

Fig.|7[ gives the a-p correlations. The origin of the 
broad bump is not fully understood. The main contri- 
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Fig. 6. The same as Fig.|^ but for '^He-p correlations. 

button results from the decay of the particle-unbound 
ground state of ^Li (J'^ = |", T ^ 1.5 MeV, Fp/T = 1, 
go = 54.3 MeV/c). Most probably, the bump contains ad- 
ditional contributions |23,EO,p7| which cannot be sepa- 
rated experimentally from the^Tl/i resonance. One contri- 
bution is expected at g ~ 16 MeV/c. It corresponds to 
the three-body decay of ^B(g.s.) — > p -I- ^Be(g.s.) — > p + 
a -|- a, where only one of the a particles is detected to- 
gether with the proton. In addition, the four-body decay 
of ioR* ^ p + 9Be(1.69 MeV) -> p -I- n -I- ^Be(g.s.) -> 
p + n -|- a + a can contribute to the broad maximum 
as was deduced at a comparably low beam energy of 40 
A-MeV 1 27 . In the present case such a contribution cannot 
be ruled out since a considerable production of 0.3 boron 
clusters per event was found in central Au+Au reactions 
at 400 ^-MeV beam energy [||. 

Due to the large Hfetime of the ^Be ground state {F — 
6.8 eV) the protons - a priori - are emitted earlier than the 
a particles. Thus, for relative momenta g ~ — 60 MeV/c 
one necessarily expects stronger final-state interaction if 
VL^a) > wl(p) and consequently R'''/R^ < 0. Indeed, the 
experimental data follow the predicted trend. However, 
since the model description does not incorporate two-stage 
decays, a reliable time difference can only be extracted for 
q values well above the resonance go - Keeping in mind that 
the sensitivity of the correlation function to a variation of 
the model parameter Ati2 decreases for increasing g, the 
deduced value is affected with a rather large error. Thus, 
the time difference derived from the a-p correlation func- 
tion will enter with negligible weight into the procedure 
used to determine the emission -tim e sequence of the dif- 
ferent particle species (cf. Sect. 1.5). 



Two resonances affect the t-d correlation function 
(Fig.H). The first one which is not resolved experimen- 
tally corresponds to the 16.75 MeV excited state of ^He 
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Fig. 7. The same as Fig. g, but for a-p correlations. The 
dark hatched area indicates the unreliable region which may 
be contaminated by doubly counted scattered particles. The 
light hatched area gives the region which is supposed to be 
strongly contaminated by secondary decays of boron isotopes. 



MeV {J"" = (|, 2; 



ond one is due to the state at 19. 

r = 2.5 MeV, g2 = 83.5 MeV/c). From the lowe"r part 
of the figure it is obvious that the apparent time delay 
is close to zero. Indeed, the model fits well the data for 
a true emission time difference (cf. table ^ which - to a 
large extent - is due to the contribution of the radial flow 
correction dTi 



The ^He-d correlations in Fig.^ appear very similar 
to the t-d correlations. The correlation function contains 
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Fig. 9. The same as Fig.|^ but for '^He-d correlations. 

one strong peak due to the 16.66 MeV state of ^Li ( J'^ = 
r = 0.3 MeV, qi = 24.8 MeV/c). Other broader 



, r ~ 5 MeV, 92 = 
(§,§)+, r = 5MeV, 



states are at 18.0 MeV (J" 
60.3 MeV/c) and at 20.0 MeV {.r - y^, 
qs = 90.2 MeV/c). Only the narrow 16.66 MeV state is 
taken into account in the model. In contrast to the t-d 
correlations, we find a clearly positive apparent time delay 
^^*He (J — 3 fm/c which increases nearly by a factor of 2 
due to the radial flow correction. 

The a-d correlation function shown in the upper panel 
of Fig.|l^ is governed by a strong resonance due to the 
narrow 2.186 MeV state of ^Li (J'^ = 3+, T = 24 keV, 
qi = 42.2 MeV/c). Another state is at 4.31 MeV {J'^ = 2+, 
r =^ 1.7 MeV, 92 = 84.2 MeV/c). Obviously, from the 
lower part of the figure one would conclude a negative 
time delay At* ^. However, for the present case the radial 



flow correction (17) overcompensates the apparent time 
difference leading to a positive value (cf. table|l|). As a 
by-product, the strong 3"*" resonance of ^Li in the corre- 
lation function can serve for an independent determina- 
tion of the q resolution. A Gaussian fit to the difference 
spectrum of true and normalized mixed yields in the re- 
gion 20 MeV/c < g <60 MeV/c delivers a dispersion of 
6vi2 = Sq/^ = 7.1 (MeV/c)//L( = 0.0057 c in good agree- 
ment with the estimate given in Sect. 2.3 . 

The a-t correlation function given in Fig. ^ is charac- 
terized by a broad maximum due to the 4.63 MeV state 
of ^Li (J'^ = I", r = 93 keV, qi = 83.4 MeV/c). Other 
states are at 6.68 MeV (J'' = f", T = 880 keV, 92 = 

116.4 MeV/c) and at 7.46 MeV (J'^ = T = 89 keV, 
93 = 126.7 MeV/c). Only the 4.63 MeV state is taken into 
account in the model. The deduced time delay At^ ^ is 
apparently negative; it is reduced by about 50 % when 
taking into account the radial flow correction (cf. table |l|). 
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Fig. 10. The same as Fig.|^, but for a-d correlations. 

The ■^He-a correlation function in Fig.|l^ is dominated 
by the 4.57 MeV state of ^Be {,r = 1^ , T = 175 keV, 
91 = 98.0 MeV/c). Here, both the deduced apparent and 
true time delays Zirgjjg ^ are found positive (table |l|) and 
differ only by about 1.5 fm/c. Similarly to the a-t corre- 
lations the small flow correction is a result of the small 
relative mass difference of the particles. 



4.3 Discussion 

Finally, the redundancy of the ten different pair-wise time 
differences obtained above allows one to fix rather reliably 
the emission-time sequence of the light charged particles. 
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Fig. 11. The same as Fig.^ but for a-t correlations. 
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Fig. 12. The same as Fig.|^, but for ''He-Q correlations. 

The optimum sequence is derived from a least-squares so- 
hition of the set of Hnear equations for the time delays 
given in table|^. (Note that, the deduced per degree 
of freedom of 1.1 does not carry quantitative information 
on the quality of the regression since the errors of the 
individual time delays are not real standard deviations.) 
The different contributions are weighted by the inverse 
squares of the given errors. If one excludes from the fit a 
few of the time differences (e.g. those which are affected 
with large errors like in the case of a-p and t-p corre- 
lations) the result changes only marginally. For a fixed 
source radius, the most probable time order of the emis- 
sion is found as follows (cf. table|2|): On the average, pro- 
tons are emitted first whereas '^He particles are emitted 
last (after about 11 fm/c). The other particles show up in 
between. Deuterons are emitted about 6 fm/c after pro- 
tons. After that, the a particles are emitted and then the 
tritons follow. However, the differences between the emis- 
sion times of these species are found in the order of 1 fm/c 
only (note: 3 fm/c lO^^s 

As a consequence of the duality of space and time coor- 
dinates, the time delays can be transformed into position 
differences. When assuming a common emission time for 
all particles, the time sequence translates into the emis- 
sion from different source radii Ri according to Eq. (^6|) 
as a consequence of the radial flow. A very similar proce- 
dure as described above for the emission times at unique 
source size leads to the average radii for the emission at 
the same time instant. The results are summarized in the 
3rd row of table^. Now, the protons, which are emitted 
earliest in the time-ordered picture, would come from the 
most extended source {Rrms = V^^o — 10 fm) whereas 
the clusters (with the exception of a particles, which ob- 
viously play a special role) are emitted from source radii 
which decrease with increasing mass. Thus, the size of the 
emission sources of the d, a, t and ■^He particles would 



be about 63%, 54%, 53% and 44% of that of the protons, 
respectively. 

This finding is in qualitative agreement with the results 
of an earlier investigation of the cluster formation pro- 
cess in central Au-|-Au collisions at 250 A-MeV. There, the 
authors compare the experimental data with predictions 
of a model which describes the hydrodynamic isotropic ex- 
pansion of an ideal nucleonic gas and the clustering by sta- 
tistical disassembly. The model predicts a breakup which, 
with elapsing time, starts at the exterior and evolves to 
smaller radii. The heavier fragments are found to arise 
from smaller source radii than the light particles. 

The OL particles do not follow completely the system- 
atic trend established above. One reason for this violation 
of the emission order with increasing mass might lie in the 
fact that this (strongly bound) particle species, at least 
partially, is made of nucleons which were initially corre- 
lated either in the target or projectile nucleus |29|] . Such 
preformed clusters would carry to some extent a memory 
of the entrance channel. 

We have studied, whether the emission time difference 
of protons and composite particles can be explained by 
calculations we have performed with the IQMD transport 
model . For a central (& < 4 fm) reaction of Ru-|-Ru we 
have investigated the distribution of the time instants at 
which the particles pass through the surface of a sphere in 
coordinate space. Only particles coming from the partic- 
ipant zone have been selected by a cut on the cm. polar 
angle | tan ©cm | > 1. Since correlation functions are sen- 
sitive to particle pairs with small relative momenta only, 
in addition, we have to demand that the velocities of both 
particles be the same. The time distributions of all the 
light charged particles are found almost symmetric and 
exhibit - within the statistical errors - identical mean val- 
ues. Alternatively, the cross check of the radial distribu- 
tions of particles with equal velocities at a unique time 
reveals no differences of the mean radii. This finding is 
not surprising since most of the transport codes generate 
composite particles by coalescence. Indeed, such a method 

Table 2. The 2nd row gives the emission times (relative to 
that of the protons) for fixed source radii Ri — Rq of d, t, 
^He, and a particles as derived from a weighted regression of 
the ten linear equations for the pair-wise time delays given in 
table |l|. The 3rd row gives the complementary information of 
different source radii (normalized to that of the proton source 
Rp = Ro) if the emission takes place at the same time instant. 





d 


a 


t 


^He 


fm/c 


6.3 ±0.8 


7.7 ±0.9 


8.5 ±0.9 


11.1 ±0.8 


R,/Rp 


0.63 ±0.04 


0.54 ±0.05 


0.53 ± 0.05 


0.44 ± 0.04 
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provides the composite particles with the same space-time 
distribution as tlie nucleons of equal velocity. 



5 Summary 

In conclusion, we have presented experimental correlation 
functions of nonidentical light charged particles produced 
in central collisions of Ru(Zr)+Ru(Zr) at 400 ^-MeV. 

For the first time an emission order of p, d, a, t, and 
■^He particles has been set up by comparing correlations of 
particles with relative momenta parallel and anti-parallel 
to the center-of-mass velocity of the pair. Collective radial 
expansion of the participant zone leads to an apparent re- 
duction of the source radius and to shifts of the emission 
times. Correcting for both effects typical time delays of a 
few fm/c were obtained. The deduced space-time differ- 
ences of the light-charged-particle emission sources allow 
two complementary interpretations. If the source radius 
is fixed the composite particles are emitted at later times 
than protons. Alternatively, if the emission time is fixed, 
the clusters are emitted from smaller sources than pro- 
tons. As a result of the duality of space and time coor- 
dinates, these two scenarios cannot be distinguished from 
each other. 
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